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We consider a general class of intermittent maps designed to be weakly chaotic, i.e., for which the 
separation of trajectories of nearby initial conditions is weaker than exponential. We show that all 
its spatio-temporal properties, hitherto regarded independently in the literature, can be represented 
by a single characteristic function <f). A universal criterion for the choice of <f) to be fulfilled by 
weakly chaotic systems is obtained via the Feigenbaum's renormalization-group approach. We find 
a general expression for the dispersion rate of initially nearby trajectories and we show that the 
instability scenario for weakly chaotic systems is more general than that originally proposed by 
Gaspard and Wang [Proc. Natl. Acad. Sci. USA 85, 4591 (1988)]. 
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The instability of deterministic motion is typically 
characterized by the growing separation of initially 
nearby trajectories, which, in turn, indicates that dy- 
namics is high sensitive to its initial state; very small 
differences at one moment in such systems can result in 
very large differences later on. For a chaotic dynamical 
system x t — T'(xo), the separation of initially nearby 
trajectories, namely Sxt = T*(:eo +Sxq) — T*(xq), evolves 
asymptotically as 



\Sx t \ ~ |&Eo|exp[A t (xo)C(*)L 



(1) 



for almost every choice of 5xq within a certain neighbor- 
hood of xq. The positive coefficient Aoo(xo) stands for 
the largest Lyapunov exponent since £(i) ~ t, and most 
of well-known chaotic systems are ruled by dispersion 
rates £(t) of this type lj-Q. After the pioneering work of 
Gaspard and Wang [4|, there has been a growing inter- 
est in recent years in understanding the so-called weakly 
chaotic (or sporadic) systems, for which the separation of 
initially nearby trajectories is weaker than exponential, 
i.e. C{t)/t — > [BI-fToj . Such systems also exhibit un- 
predictable behavior although they have zero Lyapunov 
exponent in the usual sense. 

In this paper we consider a general class of maps of 
the interval, which will hopefully be weakly chaotic, and 
we tie all of its spatio-temporal properties, including the 
map equations itself, to a single characteristic function 
<t>. The first step towards a unifying framework is based 
on determining the eigenfunctions of the Feigenbaum's 
doubling operator [Tl| 



Fg{x) = ag{g(x/a)), 



(2) 



where a is a rescaling factor. The eigenfunction that 
outlines the universality class is a fixed point of J 7 , i.e. 
Fg* = g* ■ We shall see that 



,(x) = 4>- l [<t>{x) 



(3) 



where </> _1 denotes the inverse of <j> and r is a constant. 

As far as we know, the weakly chaotic behavior 
C(t)/t — > results from the intermittent switching be- 
tween long regular phases (so-called laminar) and short 
irregular bursts [4T-flfj|]. Since the early 1980's, the idea 
that the same functional equation employed by Feigen- 
baum for studying the period-doubling cascade can also 
be used to describe intermittency and dissipative systems 
has been developed by many authors, notedly |12h16| . 
From a viewpoint more connected to the intermittency 
phenomenon, considerations done so far about universal- 
ity have been based on the scaling properties of the lam- 
inar length (l2rll5| . Our results also shed some light on 
the scaling hypothesis behind the relationship between 
laminar length and Feigenbaum's doubling operator @. 

We provide a general formula for the dispersion rate 
C(i)j one of the main results of this manuscript. But be- 
sides an attempt to outline the subexponencial instability 
in weak chaos, our goal is to provide a full description of a 
nonlinear dynamical system whenever a scaling property 
is present. Here we have two key quantities determin- 
ing the spatio-temporal properties of weakly chaotic sys- 
tems. The invariant density p(x) gives us the measure of 
concentration trajectories at each stage of intermittency, 
whereas the residence times at each ones are ruled by 
the waiting-time probability density function ip(t) of the 
laminar phase. By introducing a proper modelling of in- 
termittency mechanism, together with a renormalization- 
group approach, we establish a universal criterion for the 
choice of <f> which enables us to predict the dispersion rate 
C(t), as well as determining p{x) and ip(t). We will also 
see that these results enable us to predict the anoma- 
lous subdiffusion for spatially extended versions of such 
systems by means of a relationship between £(i) and the 
mean squared displacement cr 2 (t). 

Let us consider the general class of piecewise expanding 
maps, from [0, 1] to itself, in the form Xt+i = T(x t ) so 
that 



T{x) 



fix) 



(4) 
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with a marginal fixed point at x = 0, i.e., f(x — >• 0) = 
and f'(x 0) = 0. Although the global form of T is less 



2 



relevant here, the mechanism of intermittency also relies 
on the existence of a interval point c such that T(c_) = 1 
and T(c+) = 0, ensuring the chaotic reinjection of trajec- 
tories from (c, 1] to the laminar region near the marginal 
point. One key ingredient for the resulting weakly chaotic 
behavior is how invariant density of T behaves near the 
marginal fixed point. First of all, p(x) must be an eigen- 
function of the Perron- Frobenius operator [l7j 



v P (x) = ]T 



p[T-\x)} 
\T>[TT\ X )]\> 



(5) 



where the sum extends over all preimages TJ X of the 
point x at which the density is to be computed. Near the 
marginal fixed point i = 0we have T^ l {x) ~ x, whereas 
T~^ (x) assumes values on the interval far away from 
zero. Then one has p(x) ~ p(x)[l — f'{x)} + constant, 
resulting in the divergent invariant density 



p{x) 



1 



/'(*) 



(6) 



near x = 0, up to a positive multiplicative constant for 
non-normalizable p(x). 

The finite-time (generalized) Lyapunov exponent of 
the map T satisfying Eq. ([TJ is 



A t (x ) 



1 



J2^\T'[T k (x )}\. 



(7) 



fe=0 



In order to obtain the dispersion rate C(^)> we expect that 
the average of Lyapunov exponent ([7]) over the initial 
condition ensemble may be estimated by using a density 
function p(x, t) so that 



(Aoo) = — !— { du ( dx p(x,u)\n\T'(x)\ 
s,W Jo Jo 



(8) 



as t — > oo. To find p(x, t) we will consider the continuous- 
time stochastic model proposed in fl8j . 

^p(x,t) = ~[v(x)p(x,t)]+C(t). (9) 

The convective equation © describes the stochastic mo- 
tion of a particle initially in the laminar phase until it 
is random reinjected back to a position on the inter- 
val after reaching (crossing) the point x — c. The term 
v(x) describes the particle's velocity in the laminar phase, 
and the reinjection source term C(t) is chosen to fulfill 
conservation of quantity J Q dxp{x,t). Equation @ can 
be completely solved using the method of characteristics 
19]. Assuming uniform initial density p(x, 0) = 1 and 
considering the auxiliary functions <f>(x) and ty(x,t) as 
follow 



1 



v(x) 



dx = <p(x), ^(x,t) 



d 



(10) 



we get the general solution for the Laplace transform 
C s [p(x,t)} = p(x,s): 



p(x,s) 



1 



^{x, s) 



v{x) 1 - *(c, s) 



(11) 



The corresponding waiting-time density function tp(t) is 
given by ip{t) ~ *(c, t). Note that 4> in Eq. dTDJ) is defined 
up to an additive constant, to be suitably chosen so that 
^{c, s — > 0) = 1 due to normalization of Thus we 

set </>(c) = 1, and therefore we have (hereafter omitting 
positive multiplicative constants) 



m ~-[<j>-\t +i)]'. 



(12) 



The convective model ([9]) fulfills the ergodic properties 
of our map T noting that the average © requires that 
p(x,t) ~ Q'{t)p(x). Upon considering this asymptotic 
behavior, we must have v(x) ~ xj p{x). The dispersion 
rate £(t) can now be calculated from the general solution 
(ITT1) . and one gets 



s[l - V(s)] 



(13) 



Later we will solve Eq. (|T3|) explicitly in terms of <p. 
Besides the waiting-time density (p~2]) , the map T and its 
corresponding invariant density can also be put forward 
here in terms of 6 as follow 



T'(x) ~ l-l/x(j)'{x), p(x) ~ -x<j>'(x). 



(14) 



Now we are ready to carry out the renormalization-group 
approach by using the length of laminar motion x* , which 
is usually known for its scaling properties fl2i - [T5L [20l] . A 
detailed justification for this choice will be given here 
later on. From the definition of <j> in Eq. (p~0|) we have 
the relation r = (j){x) — </>(£*), r being the time interval 
between laminar positions x* and x. Thus the eigenfunc- 
tion of doubling operator ((2J is g*(x) =x*, as presented 
before in Eq © . The boundary conditions must be such 
that 4>'{x — y 0) = — oo, since we must have T'{x — > 0) = 1 
and positive invariant density. Then one has g*(x) ~ x, 
so that <7*(x) has the same characteristics of map T for 
x near zero, i.e. 



g»(x^0)=0, gl{x^0) = l. 



(15) 



For our first example, first notice that the eigenf unc- 
tion g*(x) of operator ([2]) can be recast in the form 
g*(x) — a<p^ 1 [(j)(x/a) — 2t}. Introducing the auxiliary 
functions h\(x) — a(f>^ 1 (2x) and h,2(x) = 4>(x/a)/2, one 
has g*(x) — hi[h,2(x) — r], admitting as possible solu- 
tions hi — (j)^ 1 and h 2 — 4>. Therefore, the characteristic 
function <fi may satisfy the relations <p{x/ a) = 2<j)(x) and 
a(f)^ 1 (2x) — 4>^ 1 {x). The first of latter two equations 
was derived in a different way in [l4j [. and it is easy to 
see that the pair admits as solution 



0(:r) 



-1/7 



2\ 



(16) 
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for 7 > 0. This result gives us the well-known class of 
Pomeau-Manneville maps, for which f(x) ~ x 1+1 / 7 near 
x = 0. From the physical point of view, the original 
Pomeau-Manneville system, i.e. 7 = 1, is paradigmatic 
since it corresponds to certain Poincare sections related 
to the Lorenz attractor [2l| . Lastly, one has p{x) ~ x -1 / 7 
and ~ (* + 1)-( 1+7 ). Our Eq. (T3J} gives us the 
corresponding weakly chaotic regimes: 



C(*) 



t(lni)" 



< 7 < 1, 
7 = 1, 



(17) 



and also ~ i for 7 > 1. 

The weakly chaotic instability Q{t)/t — > stems from 
the diverging behavior of invariant measure /z near the 
marginal fixed point, e.g. < 7 < 1 for the class 
of Pomeau-Manneville maps (Tlgl) . From Eq. ([H]) we 
have the invariant measure (J,([x, y]) ~ a;<^|^ — ducf>(u). 
Therefore, the divergence of invariant measure occurs 
provided that the characteristic function <f> obeys 



dx<f>(x) — > —00 



(18) 



as x — > 0, which is more restrictive than the boundary 
conditions (TT5)) . From Eq. (TT2"j) we also have the relation 
/ dtip(t)t ~ — J dx<p(x), and therefore Eq. (TT8l) implies 
the divergence of mean waiting time. Thus, the weakly 
chaotic behavior occurs provided that ip(t) does not de- 
crease faster than t~ 2 . In such cases, Eq. (fT3"|) can be 
solved by making use of Karamata's Abelian and Taube- 
rian theorems for the Laplace-Stieltjes transform 22|. By 
considering the general form of cumulative distribution 
function associated to ip(t), i.e., f Q duip(u) ~ 1 — l/l(i)i 7 
and l(t) being a slowly varying function at infinity, one 
obtains from Eq. (fT3"|) 



at) 



duC„ 



Ul/s) 



,1+7 



(19) 



noting that for 7 = 1 we must have l(t — > 00) = 0, thus 
excluding only the 7 = 1 border case ip(t) ~ t~ 2 (cor- 
responding to 4>{x) ~ Furthermore, noting that 
J duip(u) ~ 1 — _1 (f), one finally gets 



C(*) 



(20) 



The result (|20|) enables us to develop models with 
weakly chaotic behavior provided that the criterion (TT5)) 
is fulfilled. Thus, our second example relies on the family 
of maps for which 4>{x) behaves as 



<f>{x) ~ <fi (x) enp(x p ) 



(21) 



for /3 > 0, provided 4>o(x) does not go to zero as fast as 
cxp(— x~P) for x — > 0. This model should be understood 
as a 7 — > limiting case for the rescaling factor a = 2 7 . 
Using again the machinery we have developed, one has 



ip(t) ~ t 1 (\nt) irrespective of (f>o, resulting in 

the strong anomaly dispersion rate 



C(t) - \n 1/fi t 



(22) 



which also agrees with Eq. ([20|) . The dispersion rates 
([T7| and ([22]) are in perfect agreement with their corre- 
sponding quantities in the infinite ergodic theory [23l. l2~ij| . 
the so-called return sequences at 12 5j . Such sequences en- 
sure a suitable time- weighted average of observables that 
converge in distribution terms towards a Mittag-Lcfilcr 
distribution (25[. Note also that our Eq. © generalizes 
the invariant densities obtained in (26j for these types of 
systems. Lastly we observe that, under condition (|18l) . 
the dispersion rate (|2U|) shows that the instability sce- 
nario for weakly chaotic systems is more general than 
that originally proposed by Gaspard and Wang in 

Yet another application of the renormalization-group 
approach: weakly chaotic maps of the type ^ have been 
extensively used in the literature to model systems that 
exhibit anomalous transport, see for instance fTo| and 
references therein. The mechanism for generating deter- 
ministic subdiffusion is based on the extended version of 
map (j4]), from [0,1/2) to the entire real line, according 
to the rules f(x + N) = f(x) + N and f(-x) = -f(x), 
where N assume integer values. This results in a series 
of lattice cells with marginal points located at x — N. 
The corresponding transport properties can be under- 
stood in terms of a continuous-time random walk picture 
of this model, with probability density of waiting times 
i/j(t) near each marginal point. The mean squared dis- 
placement a 2 (t) for such model is given by [27J 



C s {a 2 (t)} 



(23) 



Since ip(s — > 0) = 1, from Eq. (fT3"j) one has 

a 2 (t)~<;(t)- (24) 

The mean squared displacements for the extended ver- 
sions of models (IT6l) and (f2"Tj) based on our results are in 
perfect agreement with those obtained respectively in [28| 
and [H, [3(| ■ Equation (|2^|) is particularly interesting be- 
cause it is a non-trivial extension for weakly chaotic sys- 
tems of a relationship typically observed in usual chaos, 
namely cr 2 (t) ~ ((t) ~ t. 

Why the laminar length does scale according to the 
Feigenbaum's doubling operator? The extension of the 
doubling operator for the class of systems discussed here 
can be understood by means of the scaling limit 



j(x) = lim a n T 2n (x/a n ), 



(25) 



since we have T n (x) ~ x + nf(x) near zero. By using 
recursion, it is simple to see that Eq. (f25j) leads to the 
doubling operator @ with g*(x) ~ x + const. /(x) and 
a = 2 7 . It is important to emphasize here that we need 
not find eigenfunctions of Eq. @ covering the whole 
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interval [0, c], but just a vicinity of zero. Thus, we can 
expand our proposal eigenfunction (J3J for r rs since 
</>(x) is singular at x = 0, resulting <7*(x) ~ a; — t/<P'(x). 
Now, from Eq. (HH), one has a;/' ~ —l/(f>', and our 
scaling hypothesis boils down simply to xf ~ /. 

We close by emphasizing that the intermittent behav- 
ior is an ubiquitous phenomenon in nonlinear science. 
Different types of intermittency (with or without exter- 
nal noise) have been studied and classified as types I-III 
intermittences pil [3l| . crisis- induced intermittency [ 32 | . 
on-off intermittency [331 . 1341 , eyelet intermittency [35l437| , 
and ring intermittency |38j. Despite their different char- 
acteristics, they all have in common the presence of some 
kind of scaling behavior. Under certain circumstances, 
some types of intermittency may even be equivalent, such 



as the type-I intermittency in the presence of noise and 
eyelet intermittency observed in the onset of phase syn- 
chronization [39j| . We believe that a major challenge 
for further works would be a general characterization of 
intermittency phenomenon having the renormalization- 
group approach as a unifying framework. 
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